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1.  INTRODUCTION 

The  purpose  of  this  study  was  to  demonstrate  that  the  surface  currents 
estimated  from  the  observed  baroclinic  structure  of  a closed  mesoscale  eddy  were 
in  fact  representative  of  the  actual  currents  encountered.  This  was  done 
following  an  experiment  in  which  the  prime  objective  was  simply  to  find  such 
an  eddy  and  determine  whether  or  not  the  structure  was  closed.  No  direct 
current  measurements  were  made  with  either  in-situ  meters  or  a towed  geomagnetic 
electrokinetograph  (GEK) . A set  of  fixes  taken  with  satellite  navigation  and 
the  dead  reckoning  (DR)  log  of  the  ship  comprised  the  only  data  available. 

To  make  matters  more  difficult  the  ship  (HMAS  KIMBLA)  was  not  fitted  with  a 
speed  log.  The  bridge  log,  hourly  sums  of  engine  revolutions  and  the 
meteorological  log  were  the  only  data  from  which  to  estimate  speed  and  DR. 

Given  these  limitations  the  analysis  is  nevertheless  worthwhile.  Currents 
were  estimated  from  the  set  of  the  ship,  i.e.  from  the  vector  differences 
between  the  calculated  progress  of  the  ship  through  the  water  and  the  observed 
progress  over  the  ground. 

2 . ESTIMATION  OF  DEAD-RECKONING  DATA 

A theoretical  model  of  ship's  speed  through  the  water  as  a function 

of  engine  revolutions  and  weather  conditions  was  constructed  to  best  fit  the 

data  from  this  particular  experiment.  Trial  19/74.  A detailed  account  of 

this  model  has  been  given  elsewhere  (Nilsson,  1977)  and  only  an  outline  need 

be  given  here.  Briefly,  the  following  effects  were  allowed  for:  retardation 

of  ship's  speed  due  to  waves  and  swell  (taking  into  account  the  relative  angle 

of  Incidence),  the  increase  in  drag  due  to  applied  rudder  (particularly  in 

beam  winds)  and  the  actual  push  received  from  winds  and  sea.  This  push 

affected  both  the  forward  speed  and  also  produced  a normal  component,  i.e. 

leeway.  The  model  was  constrained  to  use  only  those  data  normally  recorded 

by  the  ship's  officers  every  four  hours.  The  pertinent  data  consisted  of  sea 

state,  wind  (absolute)  velocity  and  swell  height  and  direction.  No  data  on 

swell  and  wave  periods  were  available.  Figure  1 shows  the  calculated  variation 

of  speed  v^  (in  the  direction  of  ship's  head)  as  a function  of  wind  velocity 

and  relative  heading.  This  function  has  been  calculated  using  the  assumptions 

that  wind,  waves  and  swell  all  come  from  the  same  direction,  that  wave  height 

h Is  proportional  to  wind  speed  w (given  by  h (m)  » 0.06  w (kn))  and  further- 
w w 
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more,  that  swell  height  is  twice  the  wave  height.  These  two  relationships 
provided  a reasonable  first  order  fit  to  the  weather  data  on  Trial  19/74 
(8-13  December  1974).  Figure  2 shows  a plot  of  the  ship's  track  for  Trial 
19/74.  A total  of  66  fixes  were  obtained  from  the  satellite  navigation  system 
during  the  course  of  the  trial  (114  hr).  The  fixes  of  significance  to  this 
work  are  lettered.  Associated  with  each  labelled  fix  is  a figure  giving  the 
number  of  hours  into  the  trial  starting  from  0001  hrs(Time  Zone  K)  on  8 T)ecember 
1974.  A list  of  the  observed  weather  data  is  given  in  the  Appendix.  Figure  3 
shows  the  dynamic  topography  of  the  area,  derived  from  Expendable  Bathythermo- 
graph (XBT)  data  taken  during  the  trial.  The  values  of  dynamic  height 
(0  - 900  m)  were  calculated  by  Andrews  (unpub.)  using  a synthetic  salinity 
profile  (Andrews  1976).  The  expected  geostrophic  currents  are  inversely 
proportional  to  the  contour  spacing. 

Now,  the  determination  of  the  model  parameters  for  ship's  speed 
rested  on  two  propositions.  Firstly,  that  the  surface  currents  over  the  legs 
CDEF  were  small  enough  for  the  net  effect  to  be  ignored.  Secondly,  that  the 
average  surface  current  over  the  leg  RS  equalled  that  over  the  leg  ST.  The 
reasonableness  of  these  propositions  can  be  judged  by  referring  to  Figure  3, 
bearing  in  mind  that  a topographic  gradient  of  10  dyn.  cm  (0/900  m)  per 
18.5  km  (10.0  n miles)  is  estimated  to  produce  a surface  current  of  0.75  m sec“^ 
- see  Section  5. 

It  was  assumed  that  the  ship's  speed  in  smooth  water,  v^,  was  linearly- 
related  to  engine  revolutions,  n,  by 


The  value  of  K was  one  of  the  most  Important  parameters  to  be  determined. 

This  was  done  as  follows:  the  weather  at  the  start  of  the  trial  was  relatively 
light  (winds  i 16  kn)  from  varying  directions  until  point  E was  reached,  after 
which  the  wind  freshened  to  30  kn  from  010°  around  point  E',  with  sea  and  swell 
reaching  1.8  m (6  ft)  and  3.7  m (12  ft)  respectively  from  the  same  direction. 
After  G the  wind  quietened  and  came  in  from  the  south.  During  legs  RST  the 

winds  were  again  less  than  16  kn  from  varying  directions.  Thus,  variations 
in  the  parameters  of  the  model  only  had  a small  effect  on  the  calculated  progress 
over  legs  RST,  but  a considerable  effect  over  the  leg  EF,  particularly  as  the 
wind  and  sea  were  almost  head  on  at  the  time. 


3 


In  view  of  the  above,  an  Iterative  procedure  was  adopted  to  fit  the 

various  parameters  to  the  data.  First  estimates  were  made  and  the  set  of  the 

ship  calculated  along  the  various  legs.  If  the  data  and  model  were  perfect, 

the  calculated  set  between  successive  fixes  would  be  due  to  surface  current 

plus  any  net  error  in  the  fixes.  This  calculated  set  has  been  termed  the 

'drift'  in  this  discussion.  In  practice,  this  will  also  contain  errors  in 

the  calculated  dead  reckoning  (DR)  data.  The  value  of  K (Eqn  (1))  was 

adjusted  to  equalize  the  magnitude  of  the  drift  over  the  legs  RS  and  ST. 

Now,  with  respect  to  the  legs  CDEF,  errors  in  calculated  ship's  speed  will 

show  up  directly  in  the  comparison  between  distance  travelled  through  the  water, 

D^,  and  distance  travelled  over  the  ground,  D assuming  no  net  effect  due  to 

surface  currents.  The  criterion  for  the  model  was  to  minimise  7(0  - D ) 

^ w g 

over  the  three  legs  CD,  DE  and  EF,  separately  and  totally.  Having  adjusted 
the  model  appropriately,  the  drift  over  RST  was  recalculated,  K readjusted  and 
so  on.  This  process  converged  after  a few  cycles  leading  to  K = 15.0  ± 0.2 
r.p.m./kn  - a value  identical  to  that  traditionally  used  by  HMAS  KIMBLA's 
officers.  The  drifts  during  legs  RS  and  ST  are  0.74m  sec“^  at  290°  and 
288°  true  respectively.  The  dynamic  topography  suggests  a mean  current  of 
0.8  - 0.9  m/sec”l  at  285°  true,  so  the  calculated  directional  data  (for  which 
there  were  no  direct  constraints)  are  most  encouraging..  However,  close 
examination  of  the  data  suggests  that  in  general  the  errors  in  calculated 
ship's  speed  is  too  large  for  useful  determination  of  current  components 
parallel  to  the  ships'  track  whenever  the  wind  speed  exceeded  20  kn.  The 
drift  component  normal  to  the  ship's  track  is  still  useful,  however,  in  cases 
such  as  that  of  leg  GH,  where  the  expected  current  is  close  to  normal  to  the 
ship's  track. 


3 . EXPECTED  ERRORS 


Before  comparing  the  observed  and  expected  current  magnitudes,  we 
’ need  to  bear  in  mind  the  expected  errors  in  the  measurements.  The  satellite 

navigation  (single  frequency)  system  should  give  fixes  with  a position  error 
' = 0.4  n.  mile.  The  larger  error  component  will  generally  be  along  the  ship's 

track,  due  to  some  uncertainty  in  ship's  speed.  However,  for  fixes  obtained 

r 

> with  the  satellite  at  high  elevation,  the  error  is  likely  to  be  greatest  in 

longitude.  The  mean  time  between  end  points  used  in  the  determination  of 
drift  values  is  3.3  hr , so  on  the  average  we  could  expect  errors 
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of  about  0.15  kn  (0 . 07 m sec~^ ) in  ea^h  component  of  drift  due  to  fixing  errors. 
Now,  with  respect  to  the  drift  component  parallel  to  the  ship's  course,  the 
error  in  calculated  speed  through  the  water  will  certainly  exceed  that,  being 
more  like  0.3  kn  (a  rough  estimate  ) in  mild  weather  and  considerably  more 
once  the  weather  rises.  A bias  of  only  1°  in  steering  also  leads  to  an 
apparent  normal  component  of  drift  0.15  kn,  so  errors  of,  say,  twice  this 
may  not  be  uncommon. 

Thus,  in  mild  weather  we  may  expect  0.2  - 0.3  kn  error  in  apparent 
current  (drift)  speed  normal  to  the  ship's  course  and  0.4  - 0.5  kn  parallel 
to  the  ship's  course.  Once  the  wind  rises,  these  errors  will  increase, 
particularly  those  parallel  to  the  ship's  course.  Under  such  conditions  it 
will  generally  be  more  appropriate  to  use  only  the  normal  drift  component. 

4.  DYNAMIC  HEIGHT  TOPOGRAPHY 

It  was  necessary  to  obtain  the  dynamic  height  topography  of  the  area 
with  respect  to  some  reference  depth  in  order  to  derive  the  expected  surface 
■ currents.  Hamon  (1965),  in  earlier  work  on  the  East  Australian  current, 

assumed  1300  m for  the  depth  of  no  motion.  This  implies  that  the  structure 
. above  that  depth  accounts  for  all  the  variability  in  surface  motion.  However, 

after  some  studies  using  free-floating  pingers  at  various  depths,  he  (Hamon,  1970) 
concluded  that  the  reference  depth  needs  to  be  greater, say  2000  m or  more. 

There  is  little  doubt,  though,  that  most  of  the  surface  currents  can  be 
adequately  described  with  reference  to  a lesser  depth.  Such  dynamic  height 
topography  has  been  obtained  in  the  past  on  a much  coarser  grid  scale  by  using 
data  from  Nansen  casts.  The  specific  volume  anomaly  is  obtained  from  water 
samples  taken  at  specific  depths  and  integrated  upwards  to  the  surface  from,  say, 
1300  m. 

> We  were  not  able  to  take  Nansen  casts  during  Trial  19/24,  but  in  any 

case  the  time  required  for  such  work  precluded  their  use  as  the  principal  tool 
V in  obtaining  a quasi-synoptic  map  of  the  dynamic  height  topography  on  as  fine 

a scale  as  we  required.  Instead,  XBT  data  were  used  to  obtain  the  dynamic 
height  structure  (0/900  m)  shown  in  Fig.  3 . The  depth  was  limited  to  900  m 
by  the  XBT  probes.  The  values  of  dynamic  height  were  obtained  in  the  following 
way:  at  the  time  the  only  XBT  probes  available  were  the  type  T4  (Plessey- 
Sippican)  nominally  limited  to  450  m depth.  Realizing  from  previous  work 


(Andrews  and  Scully-Power , 1976)  that  the  eddy  structure  extended  below  450  m, 
some  T4  probes  were  modified  by  rewinding  nearly  500  m of  the  wire  from  the 
launcher  spool  (in  each  probe)  to  the  probe  spool.  Indeed,  this  operation 
was  performed  on  ship-board  as  the  trial  progressed,  using  an  electric  drill, 
and  much  care.  Probes  so  modified  are  slightly  heavier  at  launch  and  allow- 
ance has  to  be  made  for  their  slightly  greater  rate  of  descent  (theoretically 
470  m in  71  sec  compared  to  450  m for  a normal  T4  probe).  A comparison  of 
data  taken  with  both  modified  and  normal  probes  verified  this  difference. 

A total  of  > 30  deep  probes  were  successfully  used  in  the  course  of  Trial  19/74. 
Some  of  them  gave  data  to  over  1000  m,  most  of  them  to  more  than  900  m. 

Now,  there  exists  a close  correlation  between  temperature  and 
salinity  for  subsurface  Tasman  Sea  waters  so  it  was  possible  to  calculate  the 
dynamic  height  anomalies  using  a synthetic  salinity  profile  (see  Andrews,  1976) 
from  all  the  XBT  data.  The  values  of  0(0/900  m)  from  the  modified  probes 
could  then  be  correlated  with  the  values  of  0(0/450  m)  . The  correlation  was 
extremely  close,  the  correlation  coefficient  being  0.995  from  17  profiles 
with  the  following  relationship: 

0(0/900)  = 1.334  0(0/450)  + 17.1  ± 1.2  dyn.cm  (2) 

This  relationship  was  then  used  to  transform  all  the  values  of  0(0/450)  to 
obtain  the  topography  shown  in  Fig.  3 . Furthermore,  the  very  good  correl- 
ation between  0(0/900  and  0(0/450)  indicates  that  the  values  of  0(0/450)  and 
0(0/900)  should  also  be  closely  related  to  0(0/1300). 


5.  BAROCLINIC  SURFACE  CURRENTS 


In  this  paper  we  shall  use  the  term  baroclinic  current  to  refer  to 
that  component  of  the  actual  current  which  varies  over  the  area  and  tends 
asymptotically  to  zero  (by  definition)  with  depth  - being  effectively  zero  for 
all  depths  greater  than  the  so-called  depth  of  no  motion.  The  other  component 
of  the  actual  current  is  the  barotropic  component  which  is  independent  of  depth 
and  has  to  be  vectorially added  to  the  baroclinic  current.  In  general  one 
expects  the  baroclinic  surface  currents  around  the  eddy  to  be  an  order  of 
magnitude  greater  than  the  barotropic  component. 

The  computation  of  expected  surface  currents  from  the  dynamic  height 
topography  (0/900  m)  shown  in  Fig  3 requires  an  assumption  of  the  depth  of  the 
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layer  of  no  motion.  The  dynamic  topography  has  only  been  calculated  with 
respect  of  900  m,  whereas  Boland  and  Hamon  (1970)  found  that  the  East  Australian 
Current  (EAC)  extends  to  a depth  of  at  least  2000  m.  They  adjusted  their 
baroclinic  velocity-depth  profiles  (their  Fig.  14)  to  give  the  same  velocity 
at  1200  m as  that  from  one  of  their  free-floating  pingers.  Ihis  assumes  no 
significant  barotropic  (depth-averaged)  component,  an  assumption  which  may 
be  out  by  perhaps  0.1  m sec"^ . From  their  profiles  the  actual  (baroclinic) 
surface  current  speed  can  be  related  to  the  surface  speed  with  respect  to 
900  m by 

v(actual)  = 1.15  X v(0/900  m)  ...(3) 

that  is,  approximately  87%  of  the  actual  dynamic  topography  at  the  surface 
is  due  to  the  structure  between  the  surface  and  900  m.  This  factor  has  been 
used  to  compute  the  expected  geostrophic  current  C from  the  contours  of 
Fig.  3 . 


There  is  another  consideration  with  respect  to  the  expected  surface 
current.  The  curvature  of  the  topographic  contours  adds  a centrifugal 
component  to  the  usual  geostrophic  current.  From  Von  Arx  (1962),  the  balance 
of  forces  is  given  by  (his  Eq.  4-35) 

2 

p — + grad  p - p f C = 0 ... (4) 

r h m 


where 


p is  the  density  of  sea  water, 
r is  the  radius  of  curvature  of  the  contours. 


gradj^  p is  the  horizontal  radial  component  of  the  total 
pressure  gradient  force, 

f is  the  usual  Coriolis  parameter  and 

C is  the  tangential  (meander)  current  speed, 
m 

Now  in  the  absence  of  curvature  (r  -*■“>),  this  reduces  to  the  usual  geostrophic 

equation  (Cm  C ) 

g 

Cg  = j-f  gradj^P  •••(5) 

where  C denotes  the  geostrophic  current  speed  in  the  absence  of 

s 

curvature. 
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Substituting  (5)  into  the  appropriate  quadratic  solution  to  (4), 

we  obtain 


(. 

m 


rf 

2 


(I 


....(6) 


Equation  (6)  illustrates  the  fact,  noted  by  Von  Arx,  that  in  the  presence 
of  curvature,  the  pressure  gradient  is  limited  by 


C 

g 


....(7) 


Thus 

near 

when 

by 


pressure  gradients  (and  surface  currents  in  balance)  must  bc’  small 
the  centre  of  an  anti-cyclonic  eddy,  where  r is  small.  Further, 
the  limiting  curvature  is  reached,  the  maximum  surface  current  is  given 


C 

m 


rf 

2 


....(8) 


when  the  pressure  gradient  is  such  that 
C 


g 


rf 

4 


i.e.  C can  have  up  to  twice  the  value  that  C would  have  if  there  were 
m g 

no  curvature. 


For  latitude  34°30'S,  f = 0.826  . I0~'*  ser~^  and  for  a not 
uncommon  pressure  gradient  of  1.0  dyn.  cm/km,  C = 1.21  m sic“^ . From  (7) 
the  minimum  value  the  curvature  can  have  is  given  by 


which,  for  this  example,  gives  r . =58.6  km.  Under  these  conditions,  the 

min 

actual  current  speed  would  be  C = 2.42  m sec"'.  It  is  apparent  that,  given 

m 

the  contour  spacing  and  scale  of  Fig.  3,  the  enhancement  of  the  normal  geo- 
strophic  velocity  due  to  curvature  should  be  taken  into  account  when  calcul- 
ating expected  surface  currents.  Indeed,  pressure  gradients  and  surface 
currents  around  EAC  eddies  may  often  reach  their  limiting  values.  A case 
in  point apparentl y occurred  in  March  1975.  An  intense  southward  current 
flow  offshore  from  Sydney  turned  east  at  36°  S.  This  can  be  seen  in  the 
contours  shown  by  Nilsson,  Andrews  and  Scully-Power , 1977  (their  Fig.  4a). 

The  contour  spacing  Increased  innnediately  the  stream  turned  east,  as  it  had 
to  in  order  to  maintain  the  balance  given  by  Eq.(4). 


L 


i. 


8 


' 


We  can  define 
r 


R = 


min 


where  0 < R < 1 


and  (6)  reduces  to 


C 

The  ratio 

g 


s"  ■ |o  - ^ ) 

g 


....(10) 


of  meander  to  geostrophic  current  speed  is  shown  in  Fig.  4 , 


6.  COMPARISON  OF  OBSERVED  AND  EXPECTED  SURFACE  CURRENTS 

Table  1 gives  the  drift  data  for  the  pertinent  part  of  I rial  19/74. 
The  individual  legs  can  be  related  to  the  track  of  Fig.  2 and  the  topography 
of  Fig.  3;  the  value  of  At(hr)  gives  the  elapsed  time  between  the  appropriate 
fixes.  The  'observed'  values  comprise  drift  distance  and  direction;  the 
mean  drift  speed  v^^  is  given  in  knots  and  m sec"^ . In  three  instances 
(e.g.  leg  G-H)  only  the  component  of  drift  normal  to  the  ship's  track  has 
been  used,  because  of  uncertainty  in  the  ship's  forward  speed.  These  are 
subscripted  by  . 


The  'expected'  set  has  been  calculated  by  roughly  integrating  the 
stream  flow  indicated  by  the  contours  in  Fig.  3 along  each  leg  of  the  track 
with  respect  to  time.  The  resultant  expected  drift  is  also  given  in  Table  1 . 
As  discussed  in  the  previous  section,  the  geostrophic  flow  from  0 - 900  m 
has  been  enhanced  by  a factor  of  1.15  to  arrive  at  an  expected  value  of  C . 

At  this  stage  only  the  contour  spacing  has  been  considered  - no  enhancement 
due  to  curvature  has  been  incorporated.  The  errors  in  the  individual  deter- 
minations of  observed  and  expected  speeds  preclude  using  each  value  of  vjj. 

as  a measure  of  meander  velocity  C in  order  to  calculate  individual  values 

m 

of  the  ratio  R (ref.  Eq.  (10)).  Also,  in  many  cases  the  radius  of  curvature 
varies  quite  rapidly,  so  that  only  mean  values  over  each  leg  could  be  obtained 
anyway.  It  would  be  interesting,  given  a well-defined  structure  such  as  this 
eddy,  to  correlate  a continuous  record  of  current  (such  as  might  be  obtained 
by  a GEK)  with  the  local  pressure  gradient  and  curvature.  The  residuals 
would  provide  information  on  ageostrophlc  components  and  turbulence  on  the 
scales  down  to , say,  2 km. 
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Before  comparing  values  of  drift  speed,  it  is  necessary  to  compare 
drift  directions.  The  expected  flow  directions  were  only  estimated  to 
the  nearest  s”  from  Fig.  3;  hence  the  agreement  shown  in  Table  I is  quite 
good.  The  median  value  of  the  absolute  differences  i s 6°  ( xcluding  the 
three  cases  in  which  only  the  component  normal  to  the  ship's  track  was  used). 
The  mean  difference  over  the  13  comparisons  is  + 5°,  which  is  not  significant. 
Clearly,  we  can  proceed  to  compare  the  absolute  speeds. 

For  this  experiment  the  data  must  be  treated  statistically.  The 
data  of  Table  1 have  been  divided  into  two  classes:  the  drift  values  obtained 
along  the  predominantly  straight  sections  of  the  contours  have  been  labelled 
'S',  while  those  obtained  on  the  bends  have  been  labelled  'B'.  Class  S 
contains  6 measurements  over  a total  of  27.8  hr,  class  B contains  10  measure- 
ments over  2A.6  hr.  Working  with  a sum  of  a series-of  drift  measurements 
has  the  advantage  of  reducing  the  overall  error.  For  example,  with  a 
sequential  set  of  fixes,  the  overall  fix  error  is  simply  the  relative  error 
between  the  first  and  last  fixes.  Even  if  the  series  is  not  sequential,  we 
can  expect  the  overall  error  to  quickly  become  unimportant.  A constant 
under  or  over-estimation  of  ship's  speed  with  a non-random  ship's  track  will 
continue  to  bias  the  drift  results;  but  nearly  random  errors,  such  as  we  might 
reasonably  expect  fromour  istimations  of  ship's  speed,  will  tend  to  cancel  out. 
Also,  instead  of  using  individual  values  of  v^^..  it  is  better  to  sum  the 
observed  drift  distances  within  each  class,  S or  B,  and  compare  these  sums 
with  the  sums  of  expected  drift  distances.  This  procedure  automatically 
gives  the  desired  weight  to  each  determination  obtained  over  var'ing  intervals 
of  time. 

Writing  0(d)  for  observed  drift  distance  and  E(d)  for  C^  At  for  each 
leg,  we  obtain: 

I 0(d)  - 41.0  r.mi.,  I L(d)  = 43.3  n.m.,  I At  = 27.8  hr 

S S G 

I 0(d)  = 39.2  m.mi.,  [ E(d)  = 33.0  n.mi.,  I At  = 24.6  hr 

B B B 

The  ratio  ^O(d)  / ^E(d)  = 0.95  for  class  S (the  'straight'  sections  of  the 
eddy)  and  gives  1.19  for  class  B (the  'bends').  For  class  S,  the  difference 
between  0.95  and  1.00  is  equivalent  to  an  average  over-estimation  of  ship's 
speed  of  'v  0.1  kn  (the  ship  mostly  went  around  the  eddy  with  the  current 
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stream)  or  an  over  estimation  of  the  correct  factor  relating  surface  velocity 

to  C (0-900m)  by  0.06  (i.e.  1.09  should  have  been  used  instead  of  l.lo;. 
g 

Either  error  is  quite  feasible.  The  distribution  of  the  unweighted  individual 

ratios  0(d) /E(d)  can  also  be  considered  for  each  class  of  data  and  for  set  S 

is  given  by  0.97  ± 0.2  . With  only  6 values  the  difference  from  unity  is 

obviously  not  significant.  For  class  B,  however,  the  difference  is  significant. 

Firstly,  the  result  for  class  S gives  us  some  confidence  concerning  botii  the 

estimation  of  ships'  speed  (in  mild  weather)  and  the  factor  relating  C^ftotal) 

to  C (0-900  m)  Secondly,  the  distribution  of  individual  ratios  Ofdi/Efay  tor 
g 

class  B is  given  by  1.31  ± 0.46  . With  10  observations,  tiie  mean  value 
differs  significantly  (at  5%)  from  unity.  Our  best  value  of  curvature  enhant  - 
ment  can  probably  be  obtained  by  first  reconciling  the  totals  of  .bser/cd  an<’ 
expected  drift  distances  in  class  S,  i.e.  along  the  'straights'.  We  shall  do 
this  by  reducing  the  multiplicative  factor  applied  to  C from  1.15  to  !.09. 
Whether  or  not  this  is  the  true  value  or  the  error  lies  in  the  estimat ions  of 
ship's  speed  or,  in  fact,  the  currents  are  indeed  less  than  theory  would 
indicate  by  5%,  hardly  matters  to  the  following  argument.  This  step  redtices 
^ E(d)  to  31.3  n miles.  The  ratio  '[  0(d)/^E(d)  now  gives  1.25  for  the  time- 
averaged  enhancement  due  to  curvature.  ^ 

It  can  be  shown  that  this  is  consistent  with  the  eddy  dimensions;  From 
. C 

Fig.  4,  the  ratio  1.25  for  (,  7:^^  corresponds  to  = 0.64  . The  time- 

averaged  value  of  C over  all  ® the  track  segments  listed  in  Table  I is 
8 

0.71  rn.sec"*.  This  corresponds  (from  (9))  to  ~ 

The  N - S dimension  of  the  eddy  (which  is  most  appropriate  to  the  measured 
enhancement)  across  the  maximum  current  stream  (approximately  coincident  with 
the  ship's  track)  is  about  110  km.  One  would  thus  expect  an  average  radius 
of  curvature  = 55  km,  in  excellent  agreement  with  our  value  of  ^ r^ 

One  can  consider  the  Rossby  radius  of  deformation  y as  a mesosrale 
measure  of  the  scale  of  instability  in  the  ocean.  Andrews  and  Scully-Power 
(1976)  have  shown  that 

y = k-^  (11) 

where  k is  the  horizontal  wave  number.  According  to  their  model  for  these 
eddies,  the  maximum  current  speed  occurs  at  a radius  of  0.49  x the  total  eddy 
radius,  R'  say.  Also,  k R'  is  the  first  non-zero  root  of  a first  order 
Bessel  function,  i.e.  k = 3.83/R'.  Substituting  0.49  R'  = ^r^  =54  km. 
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we  obtain  R'  = 110  km  and,  from  (11),  w=  29  km.  We  observe  that  the  value 
of  p corresponds  approximately  to  the  value  of  ^ r . . Also,  these 

calculations  are  mostly  based  on  the  smaller  north-south  eddy  dimensions, 
rather  than  the  larger  east-west  measurement.  Andrews  and  Scully-Power  used 
a slightly  larger  value  for  radius  R'  (124  km)  to  obtain  p = 32  km,  comparing 
that  with  the  value  40  km  derived  by  Godfrey  (1973  and  private  communication;. 
Indeed,  it  should  be  noted  that  the  eddy  under  study  in  this  paper  is  probably 
the  same  eddy  as  that  studied  by  Andrews  and  Scully-Power  three  months  earlier 
The  intensity  seems  to  have  diminished  by  about  15-20%  in  the  interim,  judging 
from  the  temperature  sections  (not  shown) . 


7.  THE  MEAN  BAROTROPIC  CURRENT 

The  question  of  where  the  eddy  as  a whole  is  going  is  of  considerable 
interest.  The  barotroplc  velocities  will  be  contained  in  the  observed  drift 
values  of  Table  2.  Consider  a set  of  drift  measurement  taken  all  around  the 
eddy  and  suppose  there  was  a mean  southward  barotroplc  component.  On  the 
eastern  edge  of  the  eddy  the  measured  (northward)  drift  should  be  diminished 
and  on  the  western  edge  the  drift  should  be  enhanced.  To  the  first  order, 
the  total  observed  drift  will  not  be  affected,  but  in  each  case  the  values 
of[o(d)  - E(d)J should  show  a southward  component.  Along  the  northern  and 
southern  boundaries  of  the  eddy,  the  drift  distances  would  not  change  apprec- 
iably (so  long  as  the  barotroplc  component  is  small),  but  again  there  will  be 
a southward  component. 

To  isolate  the  mean  barotroplc  component,  first  reconcile  the 

observed  and  expected  sums  of  drift  distances  in  each  class  (the  'straights, 

S and  the  'bends'  B) . As  discussed  in  the  previous  section,  this  can 

Immediately  be  done  for  class  S by  multiplying  all  the  values  of  C in  Table  1 
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by  0.95  . Then  the  values  of  C (the  expected  meander  speed)  can  be  obtained 

m 

from  each  value  of  in  class  B by  further  multiplying  by  1.25  . This  gives 

\ 0(d)  = \ E(d)  = 41.0  n mi. 

S S 

and  \ 0(d) 

B 


= \ E(d)  = 39.2  n ml. 

B 


K 
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k 

k 

I 

I 


Now  sum  the  drift  vectors  d 
in  Table  2.  We  obtain 

V ^(d)-  ¥(^) 

A-Z 


(both  observed  and  expected)  for  all  16  measures 


[ 1-6 
-.8 


n ml  for  the  north  and  east  components. 


This  net  drift  over  52.4  hr  is  equivalent  to  a mean  barotropic  component  of 
2 cm, sec”'  towards  333°  true.  We  need  to  obtain  some  idea  of  the  uncer- 
tainty of  this  measure.  The  ship's  track  from  point  J through  to  point 
S constitutes  a complete  circumnavigation  of  the  eddy  (around  the  stream  of 
rac><imum  current)  so  consider  these  8 measures  as  a separate  group. 

W obtain: 


I 0(d)  - E(d) 
J-S  " 


2.2 

3.5 


n mi 


This  net  drift  over  32.5  hr  gives  a value  for  the  mean  barotropic  current 
of  7 cm. sec”'  towards  058°  true.  Because  of  the  closed  track  and  mild  weather, 
this  is  probably  our  'best'  measure.  If  we  decrease  all  the  expected  currents 
by  10%  (a  reasonable  margin  of  error)  the  barotropic  current  is  5 cm. sec”' 
towards  074°  and  if  we  increase  them  all  by  the  same  factor  the  barotropic 
current  is  8 cm. sec”'  towards  048°.  From  this  group  of  results  we  can 
reasonably  determine  the  net  drift  of  the  whole  eddy  as  5 ± 5 cm. sec”'  towards 
045°  ± 45°.  Any  drift  west  of  north  seems  unlikely  as  the  eddy  was  already 
very  close  to,  but  separated  from,  the  continental  shelf. 


Values  for  the  apparent  drift  of  the  eddy  can  be  sought  in  another 
way:  some  areas  of  the  eddy  were  traversed  twice  during  the  trial,  as  can  be 

seen  in  Fig.  2 For  example,  XBT  data  were  obtained  in  proximity  to  points 
J and  S after  a hiatus  of  about  38  hr.  Similar  proximate  data  were  obtained 
near  points  B and  P after  a gap  of  68  hr  and  near  points  H and  N after 
30  hr.  In  no  case  did  any  difficulty  obtain  in  drawing  a quasi-synoptic  map, 
such  as  Fig  3,  Indicating  that  the  contours  themselves  were  not  appreciably  in 
motion  over  the  ground.  If  we  analyse  the  apparent  variations  of  0(0/900) 

with  latitude  over  these  sections  of  track,  we  find  that  only  the  data  near 
J and  S indicate  any  apparently  significant  meridlanal  movement  and  that 
is  small  (=  0.8  n mi,  S)  and  cannot  really  be  regarded  as  significant.  The 
apparent  movements  were: 


Area 

Time  lapse 

Movement  (N) 

B-P 

68  hr 

0 1 1 cm  sec”' 

H-N 

30 

-213 

J-S 

38 

- 1 1 1 

A. 


- 2 ± 3 
-111 
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The  uncertainties  refer  only  to  the  actual  scatter  in  the  data, 

A similar  procedure  can  be  followed  in  the  area  M-Z  to  look  for 
any  zonal  movement.  With  only  two  XBT  profiles  at  the  end  of  the  track 
(near  Z)  to  use  for  comparison,  the  result  must  again  be  regarded  with 
caution.  Both  profiles,  however,  indicate  an  apparent  shift  of  = 3.5  n mi. 
E over  41  hr,  equivalent  to  an  eastwards  movement  of  5 cm  sec“^ . 

The  data  on  the  eastern  side  of  the  eddy  were  obtained  with 
insufficient  time  lapse  (-  16  hr)  for  any  reliable  measure  and  no  satellite 
fixes  were  obtained  for  over  5 hr  immediately  preceding  point  R.  Neverthe- 
less, the  apparent  shift  is  - 5 n mi.  E.  This  leads  one  to  speculate  that 
the  eastwards  movement  may  have  been  more  pronounced  at  the  eastern  edge  of 
the  eddy  and  that  perhaps  this  was  associated  with  the  east-west  orientation 
of  the  eddy's  long  axis. 


8.  CONCLUSIONS 

Despite  limited  data  on  ship's  speed  it  has  been  possible  to  compare 
expected  currents  around  an  eddy  with  those  measured  from  the  set  of  the  ship. 
The  data  do  not  indicate  any  overall  departure  from  the  currents  expected  from 
the  dynamic  relief  estimated  to  1300  m,  once  the  theoretical  speed  enhance- 
ment due  to  centrifugal  force  is  taken  into  account.  The  time-averaged  value 
for  this  enhancement  factor  in  the  neighbourhood  of  the  eddy  curves  was 
estimated  to  be  1.25,  a value  consistent  with  the  dimensions  of  this  eddy. 

The  normal  geostrophic  current,  can  be  enhanced  by  a factor  up  to  2.0  in 
localized  areas  off  East  Australia. 

The  ship's  set  data  were  analysed  for  a mean  barotropic  component, 
estimated  as  5 ± 5 cm  sec"^,  04  5'^±  45°,  The  drift  of  the  dynamic  relief  contours 
appeared  to  be  1 ± 1 cm  sec”^  S,  5 ± 5 cm  sec”'  E with  some  indication  that 
the  eastern  front  was  expanding. 

This  type  of  analysis  holds  considerable  promise  for  the  future. 

If  the  ship  were  fitted  with  dual-axis  electromagnetic  log  or  a doppler  sonar 
in  addition  to  a satellite  navigation  system,  reliable  current  averages  over 
1 to  5 hr  periods  could  be  obtained.  Navigation  systems  rn  the  future  will  be 
able  to  fix  the  ship's  geographic  position  more  accurately  and  more  often. 

Thus  more  frequent  current  measurements  could  be  made.  If  the  ship's  position 
were  known  continuously,  a near-continuous  readout  of  measured  current  would 
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be  available.  Such  data,  unlike  those  from  a towed  GEK,  would  include  the 
barotropic  component  and  would  greatly  enhance  our  knowledge  of  oceanography. 
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THE  OBSERVED  WEATHER  DURING  TRIAL  19/74 
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Original  units  are  given  (1  m - 3.28  ft) 
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SHIP'S  SPEED.  V,  (KNOTS) 


Fig.1.  MODEL  OF  HMAS  KIMBLA'S  SPEED  (AT150RPM) 
AS  A FUNCTION  OF  WIND  VELOCrV  (KNOTS). 


152' 


FIG.  3.  DYNAMIC  TOPOGRAPHY,  D O-900m  (dyn.  cm.) 
HMAS  KIMBLA  TRIAL  19/74. 


FIG.4.  THE  RATIO  OF  MEANDER  TO 
GEOSTROPHIC  SPEED  AS  A 
FUNCTION  OF  CURVATURE  . 
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